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We explore the cosmological consequences of the recently released Union2 sample of 557 Type 
la supernovae (SNIa). Combining this latest SNIa dataset with the Cosmic microwave background 
(CMB) anisotropy data from the Wilkinson Microwave Anisotropy Probe 7 year (WMAP7) obser- 
vations and the baryon acoustic oscillation (BAO) results from the Sloan Digital Sky Survey (SDSS) 
Data Release 7 (DR7), we measure the dark energy density function f(z) = pde{z) / pde(0) as a free 
function of redshift. Two model-independent parametrization methods (the binned parametriza- 
tion and the polynomial interpolation parametrization) are used in this paper. By using the \ 2 
statistic and the Bayesian information criterion, we find that the current observational data are 
still too limited to distinguish which parametrization method is better, and a simple model has 
advantage in fitting observational data than a complicated model. Moreover, it is found that all 
these parametrizations demonstrate that the Union2 dataset is still consistent with a cosmological 
constant at la confidence level. Therefore, the Union2 dataset is different from the Constitution 
' SNIa dataset, which more favors a dynamical dark energy. 
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I. INTRODUCTION 



Since the observations of type la supernovae (SNIa) first indicated that the universe is undergoing accelerated 

■ expansion at the present stage (T], H| , dark energy (DE) has become one of the most important problems in modern 
["»>« ' cosmology. Many cosmologists believe that the cosmological constant fits the observational data well. One also has 
C*") i reason to dislike the cosmological constant since it poses the fine-tuning and cosmic coincidence puzzles [3| . A variety 

of proposals for dark energy have emerged, such as quintessence phantom fc-essence Q, tachyon Q, holographic 
. @|) agegraphic @, hessence [l(J, Chaplygin gas [ijj, Yang-Mills condensate jl2j ]. etc. 

■ A most powerful probe of DE is SNIa, which can be used as cosmological standard candles to measure directly 
the expansion history of the universe. Recently, a large sample of SNIa, the Union2 SNIa dataset [13[, was released. 
This sample consists of 557 SNIa, covers a redshift region of < z < 1.4, and is the largest SNIa sample to date. 



The Union2 dataset has been used to constrain various theoretical models [14j. However, to our best knowledge, this 
sample has not been analyzed by using a parametrization method that does not depend on any theoretical DE model, 
k^l What we shall do in this paper is just this. 

Although constraining the equation of state parameter w of DE is a popular and widely-used method to investigate 
DE, Wang and Freese [15| pointed out that the DE density pd e can be constrained more tightly than w given the same 
observational data. So in this paper, combining the latest Union2 dataset with the cosmic microwave background 
(CMB) anisotropy data from the Wilkinson Microwave Anisotropy Probe 7 years (WMAP7) observations [16| and 
the baryon acoustic oscillation (BAO) results from the Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7) (l7j . 
we measure the DE density function f(z) = Pde{z) / pde(0) as a free function of redshift. 

Two model-independent parametrization methods are used in this paper. First, we use the binned parametrization 
in which the redshifts are separated into different bins and f(z) is set as constant in each bin. Binned fits of w 
have been applied before [l8j . and similar analysis have been performed for density binning [19j . In a previous work 
[2pj . we presented a new binned parametrization method. Instead of setting the discontinuity points of redshift by 



•Electronic address: swang@mail, ustc.edu.cn| 
^Electronic address: renzhe@mail.ustc.edu.cn 
t Electronic address: mli@itp.ac.cn 



2 



hand, we treated the discontinuity points of redshift as free parameters and let them run freely in the redshift range 
covered by SNIa data. As shown in [20], this method can achieve much smaller Xmin- We shall use this new binned 
parametrization method to analyze the Union2 sample. 

Following [2lj . we also consider the polynomial interpolation parametrization in which f(z) is parameterized by 
interpolating its own values at the redshifts Z\ — i * Z max /n (i = 1, 2, ...n). Since the maximum redshift of the Union2 
sample is 1.4, in this paper we choose z max = 1.4, and consider the cases of n = 3, n = 4 and n = 5. As in [2l|, we also 
set f(z) as a constant in the range z > 1.4 where DE is only weakly constrained by the CMB data. Compared with 
the binned parametrization, the advantage of the polynomial interpolation parametrization is that the DE density 
function f(z) can be reconstructed as a continuous function in the redshift range covered by SNIa data. 

This paper is organized as follows. In Sec. II, we describe the mo del- independent parametrizations considered here 
and the method of data analysis. In Sec. Ill, we introduce the observational data and describe how they are included 
in our analysis. In Sec. IV, we present the results obtained in this paper. In the end, we give a short summary in 
Sec. V. In this work, we assume today's scale factor a — 1, so the redshift z = a^ 1 — 1; the subscript "0" always 
indicates the present value of the corresponding quantity, and the unit with c = H = 1 is used. 



II. PARAMETRIZATIONS AND METHODOLOGY 



Standard candles impose constraints on cosmological parameters through a comparison between the luminosity 
distance from observations and that from theoretical models. In a spatially flat Friedmann-Robertson- Walker (FRW) 
universe (the assumption of flatness is motivated by the inflation scenario) , the luminosity distance g?l is given by 

with 

e(z) = h(z)/h = [n r (i + Z f + n m (i + zf + (i - n r - n m )f( z )] 1/2 . (2) 

Here H (z) is the Hubble parameter, Hq is the Hubble constant, f2 m is the present fractional matter density, and tt r 
is the present fractional radiation density, given by (l6| . 

Q r = Q 7 (l + 0.227LV e// ), n 7 = 2.469 x KT 5 /^ 2 , N eff = 3.04, (3) 

where f2 7 is the present fractional photon density, h is the reduced Hubble parameter, and N e ff is the effective 
number of neutrino species. Notice that the DE density function /(z) = Pde(z) / pde(O) is a key function, because a 
DE parametrization scheme enters in f(z). 

In the following we shall parameterize f(z). Notice that the maximum redshift of the Union2 SNIa sample is 1.4, 
and in the range z > 1.4 f(z) is only weakly constrained by the CMB data. As in [2lj ]. we set f(z) to be a constant 
in the range z > 1.4, i.e., 

/(z) = /(1.4) (z>lA). (4) 

As pointed out in [2l|, fixing f(z > 1.4) can help us to avoid making assumptions about early DE that can propagate 
into artificial constraints on DE at low z. 

First, we use the binned parametrization, thus the DE density function f(z) is parameterized as, 

■^ Z ) = { a z~ x <z <z l {2<i<n) ■ ( 5 ) 

Here is a piecewise constant, and from the relation /(0) = 1 one can easily obtain e± = 1. It should be mentioned 
that there are different opinions in the literature about the optimal choice of redshift bins in constraining DE. In 
[l|j|, the authors directly set the discontinuity points by hand. In [l^, Wang argues that one should choose a 
constant Az for redshift slices. In this work, to have the maximal freedom and to have the most model-independent 
parametrizations, we do not fix the discontinuity points Zi and let them run freely in the region of < z < 1.4. As 
shown in [20|, this method can achieve much smaller Xmin- Here we consider the n = 2, 3 cases. 

Next, we use the method of polynomial interpolation to parameterize f{z). As in [2l|, we choose different redshift 
points Zi = i * z max /n (i = 1, 2, ...n), and interpolate f(z) by using its own values at these redshift points. This yields 



m = Y,f u {z ~ z tf z ~ Zt ~f~ Zt+ Vr Zn) v 

' (Zi - zi)...(Zi - Zi-i)(Zi - z l+ i)...(zi - z n ) 



i=l 
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TABLE I: Relationship between number of free model parameters n p and A\ 2 



rip 


Ax 2 (l^) 


A X 2 (2a) 


1 


1 


4 


2 


2.30 


6.17 


3 


3.53 


8.02 


4 


4.72 


9.72 


5 


5.89 


11.3 



where fi = f(zi) and z n = z max = 1.4. Based on the relation /(0) = 1, one parameter can be fixed directly (in this 
paper we fix f\ ), and only n — 1 model parameters need to be determined by the data. Here we consider the cases of 
n = 3, n = 4 and n = 5. Our parametrization is very similar to that of [2lj . 

In this work we adopt x 2 statistic [23] to estimate model parameters. For a physical quantity £ with experimentally 
measured value £ (, s , standard deviation cr^, and theoretically predicted value £ t h, x 2 is 



,2 +v,„ c,™ ^r „n ,.2 



2 _ (£obs Ctfe) /-\ 
Xf- -2 ■ 10 



The total x is the sum of all x|s, i.e. 



x 2 



= E4 w 



The best-fit model parameters are determined by minimizing the total x • Moreover, by calculating Ax = X ~Xmini 
one can determine the la and the 2c confidence level (CL) ranges of a specific model. Notice that for the la and the 
2<T CL, different n p (denoting the number of free model parameters) corresponds to different Ax 2 - Therefore, we list 
the relationship between n p and Ax 2 in table U from n p = 1 to n p = 5. 

For comparing different models, a statistical variable must be chosen. The Xmm is the simplest one, but it has 
difficulty to compare different models with different number of parameters. In this work, we will use Xmin/dof as a 
model selection criterion, where dof is the degree of freedom defined as 

dof = N- rip, (9) 

here N is the number of data. Besides, to compare different models with different number of parameters, people often 
use the Bayesian information criterion [24[ given by [25[ 

BIC = X 2 min +n v \nN. (10) 

It is clear that a model favored by the observations should give smaller Xmm/^ / an( i BIC. 



III. OBSERVATIONS 



First we start with the SNIa observations. We use the latest Union2 sample including 557 data that are given in 
terms of the distance modulus Hobs(zi) [Hj]. The theoretical distance modulus is defined as 

fhh(zi) = 51og 10 £)i(2i) + no, (11) 

where /zo = 42.38 — 51og 10 h, and in a flat universe the Hubble-free luminosity distance Dl = H^d^ (cIl denotes the 
physical luminosity distance) is 

D * ) = il + X) foWfi> (12) 
where 9 denotes the model parameters. The x 2 for the SNIa data can be calculated as 

2 /m \ " [Vobsjzj) - fith(zj;9)] 2 , . 

Xsn(.V) = 2^ Z2 ' (.i^ 

1=1 i 
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where [i bs{zi) and Oi are the observed value and the corresponding la error of distance modulus for each supernova, 
respectively. For convenient, people often analytically marginalize the nuisance parameter /i (i.e. the reduced Hubble 
constant h) when calculating x%n [IB] . 

It should be stressed that the Eq.(fT3"f only considers the statistical errors from SNIa, and ignores the systematic 
errors from SNIa. To include the effect of systematic errors into our analysis, we will follow the prescription for using 
the Union2 compilation provided in [27j. The key of this prescription is a 557 x 557 covariance matrix, Csn, which 
captures the systematic errors from SNIa (This covariance matrix with systematics can be downloaded from |27|). 
Utilizing Csn, we can calculate the following quantities 

A = ( M f 8 - ^)(Cs N h(pf S 4 ), (14) 



557 



,obs 



(15) 



557 

c= ( 16 ) 

Thus, the \ 2 for the SNIa data is [27| 

X%n=A-^. (17) 

Different from the Eq. (1X3]) , this formula includes the effect of systematic errors from SNIa. 

Then we turn to the CMB observations. Here we employ the "WMAP distance priors" given by the 7-year WMAP 
observations [l6|. This includes the "acoustic scale" I a, the "shift parameter" i?, and the redshift of the decoupling 
epoch of photons z* . The acoustic scale I a is defined as [l6| 

Za S (1 + *)^W (18) 



Here Da(z) is the proper angular diameter distance, given by 



1 f* dz' 

D A {z) = — - / -—, 19 
1 + z J E(z') 

and r s (z) is the comoving sound horizon size, given by 

, , 1 f 1/{1+z) da , . 

rJz) = -= / =, (20 

VSJo a 2 H{a)y/YTWhl^j)a 

where Of, and f2 7 are the present baryon and photon density parameters, respectively. In this paper, we adopt the 
best-fit values^ = 0.022765/i" 2 and fl 7 = 2.469 x 10 5 ^~ 2 (for T cmb = 2.725 K), given by the 7-year WMAP 
observations [Hj]. The fitting function of z* is proposed by Hu and Sugiyama (28j |: 

z* = 1048[1 + 0.00124(rj 6 /i 2 )-°- 738 ][l + gi(Q m h 2 ) 92 }, (21) 

where 

_ 0.0783(O b fe 2 )-°- 238 _ 0.560 
Sl ~ l + 39.5(fW) 0763 ' 92 ~~ 1 + 21.1(^ 2 ) 1 - 81 ' ( ' 

The shift parameter R is defined as 29] 



R(z*) = ^Q m H$(l + z*)D A {z*). (23) 
Following Pel. TtJ, we use the prescription for using the WMAP distance priors. Thus, the \ 2 f° r the CMB data is 

1 V .{.^obs _ t 
- , CMB)iJ\ x j x j 



XCMB — ( x i bS ~ x l h )(CcMB)ij( xC j bS ~ x j )' (24) 
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where Xi = (Ia,R,z*) is a vector, and (Comb)^' is the inverse covariance matrix. The 7-year WMAP observations 
fl6j give the maximum likelihood values: Ia( z *) — 302.09, R(z*) = 1.725, and z* = 1091.3. The inverse covariance 
matrix is also given in Ref. [l6| 

/ 2.305 29.698 -1.333 
( c cmb) = 29.698 6825.27 -113.180 | . (25) 
\ -1.333 -113.180 3.414 

At last we consider the BAO observations. The spherical average gives us the following effective distance measure 

m 



D v (z) = 



(l + zfD\{z) 



H(z) 



1/3 



(26) 



where Da{z) is the proper angular diameter distance given in Eq.ip]!, As in [13], we focus on a quantity 

r s {zd)/Dv (0.275), where r s is given in Eq. (f2"U|) . and Zd denotes the redshift of the drag epoch, whose fitting for- 
mula is proposed by Eisenstein and Hu [30] 

1291 (O, h 2 ) - 251 

^ = T^^^ ^ + hiShh) ^ (27) 

where 

6i = 0.313(fi m /i 2 )- - 419 [1 + 0.607(fi m /i 2 ) ' 674 ] , b 2 = 0.238(O m /i 2 ) ' 223 . (28) 

The BAO data from the spectroscopic SDSS DR7 galaxy sample galaxy sample [13] give r s (zd)/ Dy(0.275) = 0.1390± 
0.0037. Thus, the x 2 for the BAO data is, 

a _ f r s (z d )/Dy (0.275)- 0.1390 V 
Xbao-1 5^ ] ■ (29) 



IV. RESULTS 



A. Binned Parametrization 



In this subsection we will discuss the binned parametrization. As mentioned above, we will consider the cases of 
n = 2 and n = 3. For simplicity, we will call 2 bins piecewise constant pd e parametrization the ACDM2 model, and 
will call 3 bins piecewise constant pde parametrization theACDM3 model. 

First we discuss the ACDM2 model. Figure [T] shows Xmin versus redshift z for the ACDM2 model, where the Union2 
sample and the combined SNIa+CMB+BAO data are used, respectively. It is found that using Union2 alone, the 
ACDM2 model achieves its minimal Xmin ~ 529.43 when z\ — 0.162, while the best-fit value and the corresponding 
la CL of the model parameters are Q m = 0.326±£;^g and e 2 = 0.847±q Using the combined SNIa+CMB+BAO 
data, the ACDM2 model achieves its minimal Xmin ~ 530.627 when z\ — 0.158, while the best-fit value and the 
corresponding lcr CL of the model parameters are f) m = 0.279 + . I ; I2 4 and e 2 — 0.947 + . 1 ;;[qq. 

In figured we plot the evolution of f(z) along with z for the ACDM2 model. Based on the best-fit results shown in 
this figure, it is found that the Union2 dataset favors a decreasing f(z), while the combined SNIa+CMB+BAO data 
favor a slowly decreasing f(z). This means that, compared with the result given by the SNIa data alone, the result 
of the combined SNIa+CMB+BAO data is more close to the ACDM model (i.e. the cosmological constant model). 
Moreover, after taking into account the error bars, the Union2 dataset is consistent with the ACDM model at la CL. 

Next we turn to the ACDM3 model. Figure [3] shows the relationship between the Xmin an d the discontinuity 
points of redshift (z\ and z 2 ) for the ACDM3 models, where the Union2 sample and the combined SNIa+CMB+BAO 
data are used, respectively. It is found that using the Union2 dataset alone, the ACDM3 model achieves its minimal 
Xmin = 528.621 when z\ = 0.162 and z 2 — 0.552, while the best-fit value and the corresponding la CL of the model 
parameters are ft m = 0.427+°; 2 ^, e 2 = 0.665^;^ and e 3 = -0.130+^;^. Using the combined SNIa+CMB+BAO 
data, the ACDM3 model achieves its minimal x m in = 529.622 when z\ = 0.162 and z 2 = 0.421, while the best-fit value 
and the corresponding la CL of the model parameters are Vt m = 0.276t ^H, e 2 = 0.915lg and e 3 = LOgil^g^. 

In figure [31 we plot the evolution of f(z) along with z for the ACDM3 model. Based on the best-fit results shown 
in this figure, it is found that the Union2 dataset favors a decreasing f(z), while the combined SNIa+CMB+BAO 
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FIG. 1: Xrnin versus redshift z for the ACDM2 models. The left panel is plotted by using the Union2 SNIa sample alone, and 
the right panel is plotted by using the combined SNIa+CMB+BAO data. 
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FIG. 2: The evolution of f(z) along with z for the ACDM2 model. The left panel is plotted by using the Union2 SNIa sample 
alone, and the right panel is plotted by using the combined SNIa+CMB+BAO data. 




FIG. 3: The relationship between the xLn an d the discontinuity points of redshift (z± and z-i) for the ACDM3 model. The left 
panel is plotted by using the Union2 SNIa sample alone, and the right panel is plotted by using the combined SNIa+CMB+BAO 
data. The x-axis represents the redshift of the first discontinuity point z±, while the y-axis denotes the redshift of the second 
discontinuity point Z2- Notice that the light-colored region corresponds to a big an d the dark-colored region corresponds 
to a small \ 2 . Since z\ < Z2 must be satisfied, the bottom-right region of the figure is always blank. 
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FIG. 4: The evolution of f(z) along with z for the ACDM3 model. The left panel is plotted by using the Union2 SNIa sample 
alone, and the right panel is plotted by using the combined SNIa+CMB+BAO data. 
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FIG. 5: The evolution of f(z) along with z given by the Constitution dataset alone. The left panel is the case of the ACDM2 
model, and the right panel is the case of the ACDM3 model. 



data favor an oscillating f(z), which is more close to the ACDM model. Moreover, after taking into account the error 
bars, one can see that the Union2 dataset is consistent with the ACDM model at la CL. 

As a comparison, for the ACDM2 the ACDM3 model, we also plot the evolution of f{z) given by the Constitution 
dataset [32| alone in figure O As seen in this figure, for the ACDM2 model, there is a deviation from the ACDM 
model at 2a CL, and for the ACDM3 model, there is a deviation from the ACDM model at la CL. This means that 
the Constitution sample more favors a dynamical DE. Therefore, the Union2 dataset is evidently different from the 
Constitution dataset. 



B. Polynomial Interpolation Parametrization 



In this subsection we will discuss the polynomial interpolation parametrization. Compared with the binned 
parametrization, the advantage of the polynomial interpolation parametrization is that the DE density function 
f(z) can be reconstructed as a continuous function in the redshift range covered by SNIa data. As mentioned above, 
we will consider the n = 3,4,5 cases. For simplicity, we will call the polynomial interpolation parametrization of 
n = 3 the PI3 model, will call the polynomial interpolation parametrization of n = 4 the PI4 model, and will call the 
polynomial interpolation parametrization of n — 5 the PI5 model. 

First we discuss the PI3 model. Using the Union2 dataset alone, the PI3 model has a minimal Xmin = 530.102, 
while the best-fit value of the model parameters are fi TO = 0.650, — —5.063 and fs = —9.970. Using the combined 
SNIa+CMB+BAO data, the PI3 model has a minimal Xmin = 530.391, while the best-fit value of the model parameters 
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FIG. 6: The evolution of f(z) along with z for the PI3 model. The left panel is plotted by using the Union2 SNIa sample alone, 
and the right panel is plotted by using the combined SNIa+CMB+BAO data. 
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FIG. 7: The evolution of f(z) along with z for the PI4 model. The left panel is plotted by using the Union2 SNIa sample alone, 
and the right panel is plotted by using the combined SNIa+CMB+BAO data. 



are O m = 0.279, fa — 1.091 and fa = 1.434. In figureEl we plot the evolution of f{z) along with z for the PI3 model. 
Based on the best-fit results shown in this figure, it is found that the Union2 dataset favors a rapidly decreasing 
f(z), while the combined SNIa+CMB+BAO data favor an oscillating f(z), which is more close to the ACDM model. 
Moreover, after taking into account the error bars, the Union2 dataset is still consistent with the ACDM model at la 
CL. These results are similar to the results of figure 2J 

Next we turn to the PI4 model. Using the Union2 dataset alone, the PI4 model has a minimal Xmin = 530.031, 
while the best-fit value of the model parameters are f2 m = 0.368, fa = 0.485, J 3 = 0.758 and fa — 1.944. Using 
the combined SNIa+CMB+BAO data, the PI4 model has a minimal y?,„v, = 530.390, while the best-fit value of the 
model parameters are Sl m = 0.279, fa = 0.993, fa = 1.158 and fa — 1.437. In figureO we plot the evolution of f(z) 
along with z for the PI4 model. Based on the best-fit results shown in this figure, it is found that both the Union2 
dataset alone and the combined SNIa+CMB+BAO data favor an oscillating f(z), while the result of the combined 
SNIa+CMB+BAO data is more close to the ACDM model. Moreover, after taking into account the error bars, the 
Union2 dataset is still consistent with the ACDM model at la CL. 

Lately we consider the PI5 model. Using the Union2 dataset alone, the PI5 model has a minimal Xmin = 529.994, 
while the best-fit values of the model parameters are f2 m = 0.312, fa = 0.829, fa = 1.048, f± = 1.394 and f$ — 1.375. 
Using the combined SNIa+CMB+BAO data, the PI5 model has a minimal Xmin = 530.145, while the best-fit values 
of the model parameters are Q, m = 0.280, fa = 0.973, fa = 1.286, fa = 1.558 and f 5 = 0.984. In figure [HI we plot 
the evolution of f(z) along with z for the PI5 model. This figure has some subtle differences with the figure [3 but 
both the Union2 dataset alone and the combined SNIa+CMB+BAO data still favor an oscillating f(z). Again, after 
taking into account the error bars, one can see that the Union2 dataset is consistent with the ACDM model at la 
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FIG. 8: The evolution of f(z) along with z for the PI5 model. The left panel is plotted by using the Union2 SNIa sample alone, 
and the right panel is plotted by using the combined SNIa+CMB+BAO data. 



TABLE II: A comparison of parametrization methods, where the combined SNIa+CMB+BAO data are used during the analysis. 
rip is the number of free model parameters. We also include the ACDM model in this table for comparison. 



Model 


n p 




Xminf dof 


BIC 


ACDM 


1 


531.192 


0.949 


537.522 


ACDM2 


3 


530.627 


0.951 


549.616 


ACDM3 


5 


529.622 


0.953 


561.271 


PI3 


3 


530.391 


0.951 


549.380 


PI4 


4 


530.390 


0.952 


555.710 


PI5 


5 


530.145 


0.953 


561.794 



CL. Therefore, based on figure [6j figure [7] and figure Ej one can see that the polynomial interpolation parametrization 
also demonstrate that the Union2 dataset is still consistent with a cosmological constant at la CL. 

C. A Comparison of Parametrization Methods 

In this subsection we will make a comparison of the parametrization methods considered in this work. Utilizing 
the combined SNIa+CMB+BAO data, we list the x min i the Xmin/dof and the BIC for these models in table ILT1 
Based on this table, two conclusions can be obtain. First, the differences among these models' Xmin are verv small. 
Therefore, the current observational data are still too limited to distinguish which parametrization method is better. 
Second, the ACDM2 model has lower Xmin/ do f and BIC than the ACDM3 model, and the cases for the polynomial 
interpolation parametrization are similar. This hints that a simple model has advantage in fitting observational data 
than a complicated model. Moreover, we also include the ACDM model in this table for comparison. It is found 
that a cosmological constant is still most favored by the cosmological observations. This result is consistent with the 
conclusions of our paper. 

V. SUMMARY 

In this work, we explore the constraint of the recently released Union2 sample of 557 SNIa on DE. Combining 
this latest SNIa dataset with the BAO results from the SDSS measurements and the CMB anisotropy data from 
the WMAP7 observations, we measure the DE density function f(z) = Pde(z) / Pde{0) as a free function of redshift. 
To extract information directly from current observational data, two model-independent parametrization methods 
are used in this paper. By using the x 2 statistic and the Bayesian information criterion, we find that the current 
observational data are still too limited to distinguish which parametrization method is better, and a simple model has 



10 



advantage in fitting observational data than a complicated model. Moreover, it is found that all these parametrizations 
demonstrate that the Union2 dataset is still consistent with a cosmological constant at la confidence level. Therefore, 
the Union2 dataset is different from the Constitution SNIa dataset, which more favors a dynamical dark energy. 
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